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We study the lattice-spin model of RbCoBi-3 which is proposed by Shirahata and Nakamura, 
by mean field approximation. This model is an Ising spin system on a distorted triangular lat- 
tice. There are two kinds of frustrated variables, that is, the lattice and spin. We obtain a phase 
diagram of which phase boundary is drawn continuously in a whole region. Intermediate phases 
that include a partial disordered state appear. The model has the first-order phase transitions 
in addition to the second-order phase transitions. We find a three-sublattice ferrimagnetic state 
in the phase diagram. The three-sublattice ferrimagnetic state does not appear when the lattice 
is not distorted. 
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1. Introduction 

A large number of experimental and theoretical stud- 
ies have been carried out for hexagonal ABX 3 -type sub- 
stances such as CSC0CI3 and CsCoBr 3 because of the 
successive magnetic phase transitions due to the strong 
spin frustration. 1-6 ' These substances consist of quasi- 
onc-dimensional spin chains because the intrachain an- 
tiferromagnetic exchange interaction is much stronger 
than the interchain antiferromagnetic exchange interac- 
tion. Substances of this type are expressed by the stacked 
triangular antiferromagnetic Ising model. This model has 
three types of ordered phases: a partial disordered (PD) 
phase, a three-sublattice ferrimagnetic (3FR) phase and 
a two-sublattice ferrimagnetic (2FR) phase. 4 ' These or- 
dered phases are schematically shown in Figs. 1(a), 1(b) 
and 1(c). In the PD phase, the spins on one of the three 
sublattices are disordered and the spins on the other two 
sublattices are ordered antiferromagnetically. In the 3FR 
phase, the spins on the two sublattices are ordered in 
the same direction and the spins on the last sublattice 
are ordered in the opposite direction. The magnitudes 
of the sublattice magnetizations are mutually different. 
In the 2FR phase, the magnitudes of the parallel sub- 
lattice magnetizations are equal. Mekata described the 
magnetic phase transitions of these substances as a two- 
dimensional model with a nearest neighbor (NN) inter- 
action and a next nearest neighbor (NNN) ferromagnetic 
interaction, as shown in Fig. 2, by mean field approxima- 
tion (MFA). Recently, it has been found that the 3FR 
phase does not appear on these substances. 6 ' 7 ' 

RbCoBr3 is also a hexagonal ABX3-type substance. 
RbCoBr3 undergoes a structural phase transition and a 
magnetic phase transition simultaneously at 37.0 K. 8 ' 9 ' 
The structural phase transition is characterized by an 
up or down shift of the -CoBr3- chain structure along 
the c-axis. Recently, Shirahata and Nakamura have de- 
veloped a model of RbCoBr3 to explain this simulta- 
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(a) PD (b) 3FR (c) 2FR 

Fig. 1. Schematic sublattice magnetizations of (a) PD, (b) 3FR 
and (c) 2FR phases. Arrows indicate the direction and the mag- 
nitude of the sublattice magnetization. Black circle indicates that 
the sublattice magnetization is zero. 

neous structural- magnetic phase transition. 10 ' They as- 
sumed that a shift of the Co 2+ ion along the c-axis could 
be treated as a pseudo Ising spin. They introduced a 
Hamiltonian consisting of two kinds of frustrated Ising 
spin variables, that is, a pseudo Ising spin of the lat- 
tice and a magnetic Ising spin. They considered that 
the exchange interaction changed with a relative posi- 
tion between two magnetic spins. They defined the ra- 
tio of the exchange interaction to an interaction of the 
pseudo Ising spin of the lattice by K and simulated this 
model using a nonequilibrium relaxation method. They 
obtained a T-K phase diagram. The ordered phases of 
the lattice are schematically shown in Figs. 3(b), 3(c) 
and 3(d). We call these ordered phases of the lattice a 
"lattice PD" phase, a "lattice 3FR" phase and a "lat- 
tice 2FR" phase. In the lattice PD phase, the Co 2+ ions 
on one of the three sublattices are not shifted and the 
Co 2+ ions on the other two sublattices are shifted in the 
opposite direction along the c-axis. In the lattice 3FR 
phase, the Co 2+ ions on the two sublattices are shifted 
in the same direction and the Co 2+ ions on the last sub- 
lattice are shifted in the opposite direction. The magni- 
tudes of the shifts of Co 2+ ions are mutually different. 
In the lattice 2FR phase, the magnitudes of the shifts of 
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Fig. 2. Exchange interactions and sublattice of triangular lattice 
model. The NN interaction J\ and the NNN interaction J2 are 
shown. (A, B, C) denotes three different sublattices, as described 
in Appendix. 

two Co 2+ ions are equal. They found that the simultane- 
ous structural-magnetic phase transition occurred when 
the energy scale of the lattice part in the Hamiltonian 
coincided with the one of the spin part. However, there 
are two points which are not clear in their results. First, 
it seems that phase boundaries in their phase diagram 
cannot be drawn continuously near k = 0. An exam- 
ination of the intermediate phases may be insufficient. 
Second, they assumed that the phase transitions were of 
a second-order. They applied only the nonequilibrium re- 
laxation method for the second-order transition. When 
the model has the first-order transitions, we have to use 
numerical methods which are applicable in the case of 
the first-order transition such as a mixed start nonequi- 
librium relaxation method. 11 - 1 

Here, we study the lattice-spin model proposed by Shi- 
rahata and Nakamura, by MFA. Our aims are to clarify 
the above two points and to obtain a complete phase 
diagram. 

2. Model and Results 

We start with the following Hamiltonian proposed by 
Shirahata and Nakamura, 

T~L = ^lattice + ^spin, (1) 

where Wiattice is an elastic part and W sp i n is a magnetic 
part. 

^lattice is expected to take an expression of elastic en- 
ergy with regard to the lattice variables, 

^lattice = Jo(Lij ~ ^(i+l)j) 2 

i,j 

NN 

+ X! Jli^ij - L ik ) 2 , 2 ) 
NNN 

* (H) 




(a) (b) 




(c) (d) 

Fig. 3. Schematic manner of chain shifts of (a) lattice Para, (b) 
lattice PD, (c) lattice 3FR and (d) lattice 2FR phases. Black 
circles indicate the transition metal ions B 2+ . Circles indicate 
the halogen ions X~ . The X~ ions construct octahedra around 
B 2 + ions. Broken line triangles connect the positions of the B 2+ 
ions on the c-plane in the prototype structure. Here, the alkali 
ions A + are omitted for simplicity. 

Here, S(jfe> an d Y^ji)* run over NN pairs and NNN pairs 
on the c-plane, respectively. We assume for simplicity 
that the lattice variable takes +1 or —1 depending on 
whether the Co 2+ ion shifts upward or downward along 
the c-axis from the prototype structure. This treatment 
is different from the previous study where = ±1 and 
0. 10 ) However, this difference has little influence on the 
results qualitatively. The subscript i denotes a site along 
the c-axis, and j denotes a site on the c-plane. Jq , Ji 
and j\ are spring constants of the NN pairs along the 
c-axis, the NN pairs on the c-plane and the NNN pairs 
on the c-plane, respectively. Jq is positive and stronger 
than J\ and J\ because the chain structure is very hard. 
j\ is negative and J\ is positive. 

7i sp i n is defined by the following expression with the 
exchange interaction dependent on the lattice variables, 

H S pi n = - ^(jf _ <^Jo\Lij - i(i+l)j|)5ij5( i+ i)j 
NN 

~ y^y^c^i ~ — Ljk\)SjjSik ^ 

i (jk) 
NNN 

- (^2 ~ Ajf |L.y - Lu\)SijSu. 

* m 

Here, we suppose that the Ising spin variable Sij takes 
+ 1 or —1. Jq, jf and jf are the exchange interactions 
of the NN pairs along the c-axis, the NN pairs on the 
c-planc and the NNN pairs on the c-plane, respectively, 
in the prototype structure. We assume that the defor- 
mation of the lattice always decreases the magnitude of 
the exchange interaction. For simplicity, the decrease in 
the exchange interaction is supposed to be proportional 
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Fig. 4. T-k phase diagram. The broken and solid lines show that 
the phase transitions are of the first-order and of the second- 
order, respectively. 

to an absolute value of difference in the lattice variables: 
\Lij — Liiji |, and coefficients are AJp 1 2 y These param- 
eters A J^ j 2 -j control effective coupling between the lat- 
tice and the spin, jf is negative and Jf is positive. 3 -* The 
intrachain exchange interaction of RbCoBr3 is antifcrro- 
magnetic. Here, we perform the unitary transformation 
USiW = (—iySi, which changes the sign of the intra- 
chain interaction. Thus we consider that Jq and A Jq arc 
positive. Jq is stronger than jf and Jf • 

We apply to the above Hamiltonian a Scalapino-Imry- 
Pincus theory and solve it numerically. This theory ex- 
plains the experimental results of CSC0CI3 quantita- 
tively. 5 ) It is expected that this theory is valid for the 
present model. The application of this method is shown 
in Appendix. We define a ratio of J s to J L by k — J s / J L 
and use 

Jo = 5, Ji = -1, J2 = 0.1, (4a) 

J S = 5k, Jf = -k, jf = 0.1k, (4b) 

AJ ( s 0il;2) = 0.2J ( s 0;li2) . (4c) 

The obtained T-k phase diagram is shown in Fig. 4. 
When the lattice state is the lattice 2FR and the spin 
state is the PD, the phase is referred to as (lattice, 
spin)=(2FR, PD). We call the prototype structure and 
paramagnetic state the lattice Para and magnetic Para, 
respectively. The broken and solid lines show that the 
phase transitions are of the first-order and of the second- 
order, respectively. 

3. Discussion 

We obtain the T-k phase diagram of the Shirahata- 
Nakamura's lattice-spin model by MFA. The previous 
study did not distinguish the 3FR state from the 2FR 
state, and called both the 3FR and 2FR states a Ferri 
state. 10 - ) This study shows that the (Ferri, Ferri) phase 
in the previous study is a (3FR, 3FR) phase. The other 
Ferri states in the previous study are 2FR states. The 
phase boundaries in the previous study intersected at 
one point. Then, the simultaneous structural-magnetic 
phase transition from (Para, Para) to (Ferri, Ferri) oc- 
curred. However, the phase boundaries in this study does 



not intersect at one point because of the presence of inter- 
mediate (PD, Para), (Para, PD) and (PD, PD) phases. 
These phases were not observed in the previous study. In 
the case of k = 2.5 given by 

2J L - J S - AJ S , (5a) 

2J\ = Jf - Ajf, (5b) 

2J 2 L = J 2 S - AJ 2 S , (5c) 

a lattice part of the Hamiltonian coincides with a spin 
part. (See eqs.(A-4a), (A-4b), (A-5a), (A-5b), (A-5c) 
and (A-5d) in Appendix.) Then, successive simultane- 
ous structural-magnetic phase transitions (Para, Para) 
-> (PD, PD) (3FR, 3FR) occur. In the case of approx- 
imately 1.8 < k < 2.5, the phase transitions (Para, Para) 
-» (PD, Para) -» (PD, PD) -► (2FR, PD) -» (3FR, 3FR) 
occur. In the case of approximately 1.3 < n < 1.8, the 
phase transitions (Para, Para) — ► (PD, Para) — > (2FR, 
PD) — > (3FR, 3FR) occur. In the case of approximately 
k < 1.3, the phase transitions (Para, Para) — ► (PD, Para) 
-> (2FR, Para) -> (2FR, PD) -> (3FR, 3FR) occur. The 
phase transitions from (PD, Para) to (2FR, PD) and 
from (2FR, PD) to (3FR, 3FR) are also simultaneous 
structural-magnetic phase transitions. The simultaneous 
structural-magnetic phase transition occurs in a wide k 
range, although it occurs only at one point in the previ- 
ous study. That is, the simultaneous structural-magnetic 
phase transition tends to occur when the interaction of 
the lattice and that of the spin are comparable. In the 
case of k > 2.5, the behaviors of the lattice and the spin 
are reversed. 

In the case of the only lattice system, k = 0, the suc- 
cessive structural phase transitions lattice Para — ► PD 
— > 2FR occur. It seemed that the phase boundaries in 
the previous study could not be drawn continuously near 
k = 0. This is because the (PD, Para) phase did not 
appear in the previous study. The (PD, Para), (Para, 
PD) and (PD, PD) phases must appear to draw contin- 
uous phase boundaries in the whole region of the phase 
diagram. When the lattice is not distorted, the fluctua- 
tion narrows the temperature range of the magnetic PD 
phase which is obtained by MFA, but the magnetic PD 
phase does not disappear. Therefore, it seems reasonable 
to assume that the phases which include the lattice or 
magnetic PD state do not disappear, even if we consider 
the fluctuation. Because the initial configurations (Ferri, 
PD), (Ferri, Para), (PD, PD) and (Para, Para) in the 
previous study seem insufficient, the (PD, Para), (Para, 
PD) and (PD, PD) phases might disappear. It is neces- 
sary that we use additional initial configurations, such as 
the (PD, Para) and (Para, PD) phases, for the noncqui- 
librium relaxation method. 

These results show the first-order phase transitions in 
addition to the second-order phase transitions. The phase 
transitions from lattice or magnetic PD to 2FR are of the 
first-order. When the lattice is not distorted, phase tran- 
sitions from the magnetic Para to PD and from the PD 
to 2FR are of the second-order and of the first-order, re- 
spectively. 6 ^ These results are consistent. In the previous 
study, this model was examined using only the noncqui- 
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librium relaxation method for the second-order phase 
transition. It is necessary that we use methods which 
are also applicable to the first-order phase transition to 
discuss this model. 

The lowest-temperature phase is the (3FR, 3FR) 
phase. The magnetic 3FR state does not appear when 
the lattice is not distorted. 6 ) Lattice and spin degrees 
of freedom and nonequivalence interactions induce the 
appearance of (3FR, 3FR) phase. 
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Appendix: Mean Field Approximation 

We apply to eq. (1) the Scalapino-Imry-Pincus theory, 
that is, we use the following MFA, 

L A Lf ~ l B L A + l A Lf - l A l B , (Ada) 

S A Sf ~ m B S A + m A Sf - m A m B , (A-lb) 
L A L B S A Sf ~ l B m A m B L A + l A m A m B L B 

+ l A l B m B S A + l A l B m A S B (A-lc) 
-3l A l B m A m B , 
L A L A +1 S A S A +1 ~ r/ A L A L A +1 + v£ s t s (+i 



(Add) 



where, (A, B and C) denote three different sublattices, 
as shown in Fig. 2. We introduce sublattice moments as 



m 



i 



and energy along the c-axis as 



1 



(A-2a) 
(A-2b) 

(A-3a) 



N 
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where, Nm is the number of spins in the chains, and A 
denotes A, B and C. The Hamiltonian for the A sublat- 
tice is written in the above MFA and can be divided into 
the lattice part and the spin part as, 

W&A - E L t L ti + £ L A , (A4a) 



<fa = J S E + H* A E St, (A4b) 



where, 



J lA =-(2J L + AJ s r^), (A-5a) 
H lA = - {6J 1 L (Z B + Z c ) + 12J 2 L ; A 

+ 3AJ 1 s (Z B m B + l c m c )m A (A-5b) 

+ 6AJl(m A ) 2 l A } 1 

J sA = - (J S - AJ S + AJ$nt), (A-5c) 
J ffS A = _ { 3 ( jS_ AJ s )(m B + m c ) 

+ 6(J 2 S - A J 2 s )m A 
+ 3Aj!(l B m B + l c m c )l A 
+ 6AJl(l A fm A }. 



(A-5d) 



Equations (A-4a) and (A-4b) can be solved by the well- 
known exact solution of a one-dimensional Ising model. 
We obtain the following self-consistent equations, 



l A = 



sinh(/3i? 1 



m 



Y / exp(-4/3J LA ) + sinh 2 (^ff LA ) ' 

smh(f3H sA ) 
v /exp(-4/3J sA ) + sinh 2 (/37J SA ) ' 



(A-6a) 
(A-6b) 



and 



tanh(2/3J LA ) 



cosh([3H l 



sinh(2/3J LA )A/exp(4/3J LA )sinh 2 (/3i? LA ) + 1 



(A-7a) 



n A - 
Vs - 



tanh(2/3J sA ) 



cosh f3H sA 



sinh(2/3 J sA ) ^/exp(4/3J sA ) sinh 2 (/3ff sA ) + 1 

(A-7b) 

where, [3 = k^T . 

The partition functions of the lattice chain and spin 
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chain are 
Z L * =cxp(/?J LA )cosh(/3ff LA ) 

+ yj exp (2/3 J LA ) sinh 2 ((3H LA ) + exp ( - 2(3 J hA ) , 

(A-8a) 

Z sA =exp(/3J sA )cosh(/3 J ff sA ) 

+ v /cxp(2/3J sA ) sinh 2 (/3if SA ) + exp(-2/3J sA ), 

(A-8b) 

respectively. The free energy F is given by 

F = -N f31n(Z hA Z hB Z L ° Z sA Z sB Z s ° ) + E int , (A-9) 
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where, N is the numbers of three chains and 
E int =| {AjS^ + V B V B + ^ s c ) 

+ 12 Jl(l A l B + l B l c + l c l A ) 

+ 12 J 2 L [(Z A ) 2 + (Z B ) 2 + (Z c ) 2 ] 

+ 6(jf - Ajf )(to a to b + m B m c + m c m A ) 

+ 6(J| - AJ|) [(m A ) 2 + (m B ) 2 + (m c ) 2 ] 

+ 18 A jS(l A m A l B m B + l B m B l c m c + l c m c l A m A ) 

+ 18AJ 2 S [(; A ) 2 (m A ) 2 + {l B f(m B f + (l c ) 2 (m c ) 2 ]}. 

(A-10) 

The i?i n t is necessary to correct the dual computation of 
the interaction. 



